Gravitational Lensing Analyzed by Graded Refractive Index of Vacuum 
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We found strong similarities between the gravitational lensing and the conventional optical lensing. The 
similarities imply a graded refractive index description of the light deflection in gravitational field. We got 
a general approach to this refractive index in a static spherically symmetric gravitational field and obtained its 
exterior and interior solutions exactly through the general relativity. In weak field case, the two solutions come to 
a simple unified exponential function of the gravitational potential. With these results, the gravitational lensing 
can be analyzed in a convenient optical way. Especially, the long puzzling problem of the central image missing 
can be solved easily. We also pointed out that the graded refraction property of the gravitational spacetime is 
related to the vacuum influenced by the gravitational matter. 
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I. INTRODUCTION 

Gravitational lensing is an effect predicted by general rel- 
ativity [1, 2, 3]. It has been now a powerful tool for the 
study of the so intractable problems in astrophysics and cos- 
mology such as the value of Hubble constant f^, the physics 
of quasars |5], the mass and mass distribution of a galaxy or 
agalaxy cluster[2], the large scale structure of the universe 
EH], the existence of dark matter i [H, the nature of 
dark energy lfTll[l2l] and so on. 

It is told that, the light deflection of a gravitational lens is 
caused by the gravity or in Einstein's words the curved space- 
time, while the light deflection of an optical lens is caused by 
the variance of refractive index of the medium. But is there 
any deeper similarity between the two lenses besides the de- 
flection of light? Can we analyze the gravitational lensing as 
the conventional optical one? 

In fact, people have long tried an optical description of the 
lensing effect. In 1920, Eddington 1 13] suggested that the light 
deflection in solar gravitational field can be conceived as a 
refraction effect of the space in a flat spacetime. The idea 
was further studied by Wilson Il4l . Dicke OJ], Fehce ifiell . 
and Nandi et al. 

Recently, fliis thought of fight 
deflection has been investigated further by Puthoff [20, 211], 
Vlokh [22], Ye and Lin [23] etc. In Puthoff's paper, the light 
deflection is related to the vacuum polarization in gravitational 
field. Vlokh discussed further this light refraction effect. Ye 
and Lin used the refractive index to simulate the gravitational 
lensing. 

The purpose of this paper is, on the basis of the general rela- 
tivity, to find the similarities between the gravitational lensing 
and the optical lensing, and to analyze the former in a sim- 
ple optical way. The paper is organized as follows: in Sec. II, 
we point out the strong similarities between the two lenses and 
suggest a refractive index analysis of the gravitational lensing; 
in Sec. Ill, we find the exact solutions of this refractive index 
for outside and inside the lens matter system respectively; in 



Sec. IV, we make a weak field approximation; in Sec. V, we 
apply the obtained result to the problem of gravitational lens- 
ing, especially to the central image missing problem puzzled 
physicists for a long time ll24ll ; in Sec. VI, we make a discus- 
sion on the refraction property of the gravitational spacetime; 
finally in Sec. VII, we draw our conclusions. 



II. SIMILARITIES BETWEEN THE TWO LENSES 
A. Similarity in Fermat's principle 

Landau and Lifshitz have derived from the general relativity 
the Fermat's principle for the propagation of light in a static 
gravitational field as follows L25i] : 



'J 800 



'l^dl = 0, 



(1) 



where dl is the length element of the passing light measured 
by the local observer, goo is a component of the metric tensor 

— 1/2 

•^oo corresponds to an element of optical path length, 
goo - dt/dT, where dr represents the time interval mea- 
sured by the local observer for a light ray passing through the 
length dl, while dt is the corresponding time measured by the 
observer at infinity. Eq. (1) could then be rewritten as 



/dt dl 
ds = 0, 
dr ds 



(2) 



where ds is the length element measured by the observer at 
infinity, corresponding to the local length dl. 

For an optical lens, the light propagation satisfies the con- 
ventional Fermat's principle 



6 I nds = 0, 



(3) 
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where n is the refractive index of the medium. 

The similarity between Eq. (2) and Eq. (3) indicates that 
if we set the scale of length and time at infinity as a stan- 
dard scale for the whole gravitational space and time, the light 
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FIG. 1: Light deflection in a gravitational lens. 



propagating in a gravitational lens could then be regarded as 
that in an optical medium with the refractive index being 



dt dl 
dr ds' 



(4) 



where dt/dr relates to the curved time and dl/ds relates to the 
curved space. 



B. Similarity in light deflection formula 



Since 



rde 

tan/? , 

dr 



associating with Eq. (8) reaches 



d6: 



dr 



y \ no»'o I 



(9) 



(10) 



Eq. (10) and Eq. (6) show another strong similarity between 
the two lenses, which once more indicates that there could be 
a refractive index analysis of the gravitational lensing. The 
similarity tells us that this refractive index can be figured out 
through the following two equations: 



dR 



dr 



R 



r^/B(R) 



(11) 

(12) 



For the detailed derivation of the above two equations, see 
the APPENDIX. 



First we consider the light deflection in a gravitational lens. 
For a static and spherically symmetric lens matter system, the 
metric has the standard form 



dT^ = B{R)c^df- - A(R)dR^ - R^(d6^ + sin^Odcf)^), 



(5) 



where R is the radial coordinate of the metric. 

The light deflection is shown in Fig. 1, where curve PqP rep- 
resents the light ray, M is the mass of the lens matter, is the 
angular displacement of the coordinate radius R, j3 is the angle 
between the coordinate radius R and the tangent of the ray at 
point P. The general relativity gives the angular displacement 
as follows i26il27il2i: 



d0-- 



dR 



Rl^lMR) 



(6) 



R/^fWRj 



1 



where Rq represents the radial coordinate at the nearest point 
Po. 

Now we consider the light deflection in a medium of a 
spherically symmetric refractive index n. According to the 
Fermat's principle, the light propagation in such a lens satis- 
fies the following relation ||29(] : 



nr sin yS = constant, 



(7) 



where r, dififering from the above said coordinate radius R, 
is the distance from the light to the center of the lens. The 
relation can be rewritten as 



nr sm/3 — noro. 



(8) 



where ro and «o represent the radial distance and the refractive 
index at the point closest to the center respectively. 



III. EXACT SOLUTIONS OF THE REFRACTIVE INDEX 

A. Exterior solution 

The coefficients A(R) and B(R) in Eqs. (11) and (12) 
can be obtained from the Schwarzschild solutions. The 
Schwarzschild exterior solution (R > Rl, Ri is the radial co- 
ordinate at the surface of the lens matter system) gives ll28ll : 



2GM 



B(R) = 1 - 



Rc^ ' 



(13) 
(14) 



where G is the gravitational constant, c is the velocity of light 
in vacuum without the influence of gravitational field. 
Substituting Eq. (13) into Eq. (11) gives 



2GM 
Rc2 



(15) 



where is a constant. Since in the Schwarzschild metric, the 
ratio r/R 1 at infinity, we have 



CM 



Thus Eq. (15) gives 



R 

'^2 



2GM 

Rc^ 



+ 1 



cm] 



or 



R 



(16) 



(17) 



(18) 
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Then through Eqs. (12), (14) and (17), we get the exterior 
refractive index 



1 



1 / /i 2GM , 1 Cm\ /i 2GM 
2\y^ Rc^- ^ Rc- j yj Rc^ 

or through Eqs. (12), (14) and (18), we get 



(19) 



(20) 



which is exactly in agreement with that given by Felice in a 
different way llql . 



B. Interior solution 

The Schwarzschild interior solution (R < R^) gives ll28ll 

-1 



Jr 



2G 

R 



B{R) = exp 
AnR^piR) 



(21) 



M(R) 



1 - 



2GM(R) 



Rc^ 



dR\i22) 



rR 1 

where M{R) - 4nR~p{R)dR, p{R) is the mass density, and 
for an ordinary lens matter system, the pressure p{R) = 0. So 



r'- 2GM{R) 
2GM(Rl) 



2GM(R) 



exp 



•JRl 



R^c^ 



1 - 



Rc^ 
2GM{Rl) 



Rc^ 



dR 



dR\. (23) 



Substituting Eq. (21) into Eq. (11), and considering the re- 
lation between the radial coordinates ri and Rl at the surface 
determined by Eq. (17), we obtain 



Ri 



^ _ 2GM{Rl) ^ ^ GM(Rl)] 



Rlc'- 



Rlc^ 



exp 



rRi 
^« R^\ 



dR 



ICMjR) 
Rc^ J 



(24) 



Then combining Eqs. (12), (23) and (24) gives the interior 
refractive index 



R 

n - ni — exp 
Rl 



•^'^ R^l 



dR 



exp 



f r-Ri 
Jr 



GM(R) 



1 - 



2GM(R) 



Rc^ 



2CM(R) 
Rc- 

-1 ^ 
dR 



where 



/ f, 2GM(J 



2CM(Rl) . 1 CM(Rl} 
^ — + i — — r — T — 

^ Rlc- 



2GM(Rl) 



(25) 



(26) 



is the refractive index at the surface. 



IV. WEAK FIELD APPROXIMATION 

For an ordinary lens matter system, the gravitational field is 
not extremely strong, i.e., GM/Rc^ or GM/rc^ « 1, then we 
have: 

The exterior refractive index 



2GM\ 



n s exp I — — I = exp 



-2Pr 



(27) 



where P,- - -GM/r is the gravitational potential at position r 
outside the lens matter system. 

The above expression of graded refractive index has been 
verified in the problem of light deflection in the solar gravita- 
tional field 112311 . with the result being in agreement with that 
given by the general relativity |3, 28J and the actual measure- 
ments 1 30]. 

And the interior refractive index 



n = exp 



2GM(ri,) 



r 



2GM(r) 



dr 



exp 



-2Pr] 



(28) 



where Pr - -GM(ri)/ri - J^''' [GM(r)/r^]dr is the gravita- 
tional potential at position r inside the lens matter system. 

In general, we have the refractive index profile both outside 
and inside the lens matter system in weak field case as follows: 



n - exp 



-2Pr 



(29) 



The above result is derived from a single static spherically 
symmetric lens matter system. For a multi-body system, the 
total gravitational potential will be the superposition of each 
potential; therefore, the refractive index can be expressed as 



n = 



-2Pr 



= exp 



-2{Prl+P,.2+Pri + ---) 



n I nana ■ ■ ■ 



(30) 



where P^i , Pi-2,Pi3, ■ ■ ■ and ni , n2, "3, ■ ■ ■ are the gravitational 
potential and the corresponding refractive index caused by 
each gravitational body respectively. This expression may be 
extended to arbitrary distributed matter systems. Fig. 2 shows 
such an example, where the brighter cyanine represents the 
higher value of refractive index, and the closed curves are the 
isolines — the denser the lines, the quicker the change of re- 
fractive index. 



V. APPLICATION TO GRAVITATIONAL LENSING 

Eq. (29) and Eq. (30) are the main results of this paper, 
which provide a convenient optical way to describe the effect 
of gravitational lensing. Considering a source S and a lens L 
of mass M, the light emitted from S is bent due to the gravi- 
tational field of the lens. The bent light could be figured out 



FIG. 2: Refractive index profile of a gravitational lens composed of 
three celestial bodies of different mass. 
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FIG. 4: Ray tracing results for the central imaging, (a) Tracing the 
central imaging ray for a given lensing system; (b) Tracing the central 
imaging rays for lenses of difi"erent mass. 
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FIG. 3: A computer simulation of gravitational lensing. 



through Eq. (7) and Eq. (29). Drawing the extension line of 
the light from the observer O, the apparent (observed) posi- 
tion of the source image could then be found out. The result 
is shown in Fig. 3, where /i, h represent the upper and lower 
images respectively. 

The refractive index given by Eq. (29) could also be applied 
in studying the formation of the central image, which is pre- 
dicted by the general relativity but not observed in almost all 
known cases of gravitational lensing. This problem has puz- 
zled people for many years ll24ll . 

For the central imaging, the refractive index profile inside 
the lens matter system should also be considered. As a model 
for discussion, we suppose a matter system (for example, a 
galaxy or a cluster of galaxies) of radius ri with a density 
distribution given by 



P=Pc 



1-1^ 



(31) 



where pc is the central density of the system, < r < r^, 
^ > 0. The density p decreases with the distance r from the 
center of mass; the decreasing rate depends on the parameter 
k. This model gives the distribution of gravitational potential 



Pri - -4kPcG 




6 {2 + k){3+k)\rL 



(32) 



'•'H33) 



for outside (r > ri) and inside (r < ri) the matter system 
respectively. 

The refractive index profile outside and inside the lens mat- 
ter system then reads 



«o = exp 
n, = exp 



SnpcG 



2(2 + k) 



&■ 3(3 + k) r 



1 1 
6 " (2 + k)(3 + k) Xn 



















II 









(34) 



.(35) 



Fig. 4 (a) shows a ray tracing result for the imaging of a 
gravitational lens with the above described refractive index 
profile. In the figure, only three paths of ray (indicated by the 
three thick lines) could pass through the observer (9, forming 
the upper, lower and central images respectively. From the fig- 
ure, we find that, the larger the distance OL from the observer 
to the lens body, the closer the central imaging ray to the lens 
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center. In the same way, the larger the distance SL from the 
source to the lens body, the closer the central imaging ray to 
the lens center This can be easily seen by interchanging the 
source and observer in the figure. 

Fig. 4 (b) shows the dependence of the ray position of the 
central image on the lens mass. The lens mass is expressed 
by M = jKrjjDck/(3 + k), which can be derived from Eq. (31). 
The four curved lines in the figure represent respectively the 
four central imaging rays in four lenses of different mass. The 
mass ratio of the corresponding lenses is 2 : 3 : 4 : 5 from 
bottom to top. We find that, when the mass M increases, the 
central imaging ray will be closer to the lens center. 

This imaging characteristic of Fig. 4 can be used to explain 
the missing of the central image of gravitational lensing. In 
practical observations, the distances OL, S L and the mass M 
are all in astronomical scale; therefore, the light of central 
imaging is extremely close to the lens center. For a lens mat- 
ter system denser in the center, the possibility of the central 
imaging light to be blocked increases. The relatively longer 
inner path of the central imaging ray adds the possibility of 
scattering and absorption on the way. Besides, the central im- 
age, if still exists, is relatively faint compared with the bright 
lens core. All these factors lead to little chance of finding the 
central image. This accounts for the central image missing in 
almost all observed cases of gravitational lensing. 



VI. DISCUSSION ON THE REFRACTION PROPERTY OF 
THE GRAVITATIONAL SPACETIME 

Though we have successfully given a refractive index de- 
scription of the gravitational lensing and solved a puzzling 
problem of it, there is still a doubt about the refraction prop- 
erty of the gravitational spacetime as shown in Eq. (4). The 
doubt can be eliminated if we consider that, around the grav- 
itational matter, there exists a special medium — vacuum, 
which may have a graded refractive index in gravitational 
field. 

A graded vacuum refractive index needs at least an influ- 
ence of matter or field on the vacuum. It is exhilarating to 
see that the recent theoretical and experimental progresses in 
quantum vacuum have provided a strong support to such an in- 
fluence. Ahmadi and Nouri-Zonoz pointed out that under the 
influence of electromagnetic field, vacuum can be polarized, 
which has led to astonishingly precise agreement between pre- 
dicted and observed values of the electron magnetic moment 
and Lamb shift, and may have effects on the photon propa- 
gation ifJUl . Rikken and Rizzo considered the anisotropy of 
the optical properties of the vacuum when a static magnetic 
field Bo and a static electric field Eq are simultaneously ap- 
plied perpendicular to the direction of light propagation Il32ll . 
They predicted that magnetoelectric birefringence will occur 
in vacuum under such conditions. They also demonstrated 
that the propagation of light in vacuum becomes anisotropic 
with the anisotropy in the refractive index being proportional 
toBoxE(). Dupayseffl/. studied the propagation of light in the 
neighborhood of magnetized neutron stars. They pointed out 
that the light emitted by background astronomical objects will 



be deviated due to the optical properties of quantum vacuum 
in the presence of a magnetic field |33]. 

The fact that the light propagation in vacuum can be mod- 
ified by applying electromagnetic fields to the vacuum indi- 
cates that vacuum is actually a special kind of optical medium 
ll3Tl[33ll . The similarity between the vacuum and the dielectric 
medium implies that vacuum must also have its inner struc- 
ture, which could be influenced by matter or fields. Actu- 
ally, the structure of quantum vacuum has been investigated 
in quite a number of papers recently ll34l l35i l36ll . 

Besides the electromagnetic field, the existence of matter 
can also influence the vacuum. For example, the vacuum in- 
side a microcavity is modified due to the existence of the cav- 
ity mirrors, which will alter the zero-point energy inside the 
cavity and cause an attractive force between the two mirrors 
known as Casimir effect ll37i[38i] . which has been verified ex- 
perimentally (HI, S . 

The above said makes it natural to suppose that the refrac- 
tive index of vacuum can be influenced by the gravitational 
matter This thought is also supported by the theory of quan- 
tum fields, which tells us that particles are actually the excited 
vacuum. And it is unthinkable that the excited vacuum will 
exert no influence on the vacuum around. 

In brief, we could relate the graded refraction property of 
the gravitational spacetime to the vacuum influenced by the 
gravitational matter. 



VII. CONCLUSIONS 

We have shown two strong similarities between the gravi- 
tational lensing and the conventional optical lensing: one is in 
the Fermat's principle, the other is in light deflection formula. 
The similarities indicate a graded refractive index analysis of 
the gravitational lensing. We derived through the general rel- 
ativity and the Fermat's principle the general expression of 
this refractive index for a static spherically symmetric grav- 
itational field. From this expression and the Schwarzschild 
exterior and interior solutions, the exact refractive index pro- 
file is obtained. In weak field case, we got a simple unified 
exponential function of the gravitational potential for the cal- 
culation of the refractive index profile outside and inside the 
lens matter system. Since the derivation is based on the gen- 
eral relativity, there is no inconsonance between the general 
relativity and our result. By using the obtained result, we in- 
vestigated the gravitational lensing in a conventional optical 
way. Some results from computer simulations interpreted the 
long puzzling problem of the central image missing in a clear 
way. We also suggested the graded refractive index of the 
gravitational spacetime be interpreted as that of the gravita- 
tional vacuum. We hope our work will be a useful means of 
gravitational lensing and be a positive stimulus to the vacuum- 
based investigation of gravitation. 
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APPENDIX: DERIVATION OF EQS. (11) AND (12) 



Suppose the relation between R and r is 
dR _ dr 



(A.l) 



where f{R) is a function of R. 

Now we consider the light path Pco^'o in Fig- 5, where Pq is 
the closest point to the gravitational center Combining Eqs. 
(6), (10) and (A.l), we get 



f(R) 



R/^/B{R) 



Rq/ Vbw J 



(A.2) 



At the infinite point P^o, the space-time becomes flat. So 
we have 



^JA(R^) = 1, 


(A.3) 


^JB(R^) = 1, 


(A.4) 






^1 -1 

J loo — -1, 

dr 


(A.5) 


Roc J 

r CO 


(A.6) 


Woo = 1. 


(A.7) 



According to Eqs. (A.l), (A.3), (A.5) and (A.6), we have 

f(Roc) = Rco/ ^JA(Rco). (A.8) 
Applying Eq. (A.2) to the infinite point Poo reads 



fiRo.) 



Roo/^^B(R^ 



- 1 = 



- 1. 



/?o/VW) 

(A.9) 

Substituting Eqs. (A.4), (A.6), (A.7), (A.8) into the above 
equation gives 



"on) 



R,,l y/B(Ro). 



(A. 10) 



Next we consider another light path P^Pq in the gravita- 
tional field. Now the closest point to the gravitational center 
is not Po but P^^ , then Eq. (A.2) becomes 



R/yfMR) 
f(R) 



R/^m 



\[K/ylB(Ro) 
Similar to Eq. (A. 10), we can get 



'o'o 



1. (A.ll) 



(A.12) 




FIG. 5: Two light paths in a gravitational field. 
Applying Eq. (A. 1 1) to the point Po reads 
Ro/^/A(Ro) 



f(Ro) 



- 1 = 



Combining Eqs. (A.IO), (A.12) and (A.13) gives 



/(Po) = Ro/ y/MRo). 



(A.13) 



(A. 14) 



In the same way, we can obtain the following relation for 
every point P in the gravitational field: 



/(P)=P/Va(P)- 



(A.15) 



Substituting the above equation into Eq. (A.l) gives Eq. 

(11) . Combining Eqs. (6), (10), (A.IO) and (11) gives Eq. 

(12) . 
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